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Abstract. Highly collimated jets found in AGB and post-AGB stars are expected to play an 
important role for shaping planetary nebulae. Recent VLBI observations of H2O maser sources 
have revealed that some of the spatio-kinematical structures of H2O maser sources exhibit stellar 
jets with extremely spatially and kinematically high collimation. Such stellar H2O maser jets, 
so-called "water fountain" sources, have been identified in about 10 sources to date. Here we 
review recent VLBI observations and relevant observational results of the water fountain sources. 
They have revealed a typical dynamical age and the detailed kinematical structures of the water 
fountains, possibility of the existence of "equatorial flows", and the evolutionary status of the 
host stars. The location and kinematics of one of the water fountain sources in the Galaxy is 
also revealed and shown here. 
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1. Introduction 



A stellar jet appearing at the final stage of stellar evolution is a manifestation of the 
important phenomena to form asymmetric planetary nebulae (PNe). Based on the large 
~H ■ incidence of bipolar and multipolar optical morphologies and the common presence of 
point-symmetry found in a sample of young PNe, [Sahai fc Trauger (1998) concluded that 



collimated jets have been operational in a large fraction of PNe, and expected to be one 
, of the major factors shaping these objects. Similar morphologies have been found in a 
t-H ' survey of pre-PNe (jSahai 2004p . leading to the current paradigm that such a collimated 
\ jet had already been launched at an earlier phase of stellar evolution - the asymptotic 
giant branch (AGB) phase - and is frequently observable in molecular line emission. 
f**«« I Subsequently an optically visible PN is seen at the post-AGB phase. 

There are candidate sources that show us the earliest stage of stellar jet emergence. 
, "Water fountains" are the most promising candidates; they exhibit extremely-high- velocity 
flows traced by H2O maser emission. The outflow velocity is much faster than the typ- 
ical expansion velocity of circumstellar envelopes of AGB stars found in 1612 MHz OH 
masers (10-25 km s _1 , e.g.. Ite Lm tcl Hckk ert et al. 1989|) . The H2O maser emission ex- 
hibits high collimation in its both morphology and kinematics. This suggests that high 
density gas is supplied from the vicinity of the stellar surface to the jet tips where H2O 
maser emission is excited (e.g., [Elitzur et al. 1989|) . Because VLBI observations have been 
made for some of these water fountains, it is possible to undertake statistical analysis 
to estimate typical values of their physical parameters, lead to understanding the final 
evolution of the circumstellar envelopes and the formation mechanism of stellar jets. 

Note that high precision astrometry of H2O masers in the water fountains provides 
great opportunity to directly estimate the locations and kinematics of the water fountains 





Stellar molecular jets traced by maser emission 
Table 1. Parameters of "water fountain" sources. 



1 







Length 


Al/i os t 


tjctjf 


Distance 




IRAS name 


Other name 


(arcsec) 


(km s 


) (year) 


(kpc) 


Reference^ 


15445—5449 


OH 326.5-0.4 




91 






4 


15544—5332 


OH 325 8—0 3 




40 






4 


16342—3814 


OH 344.1+5.8 


2.4 


240 


~100 


2 


3, 15, 16, 17, 18 


16552-3050 


GLMP 498 




170 






19 


18043-2116 


OH 0.9-0.4 




204 






4 


18139-1816 


OH 12.8-0.9 


0.11 


55 


70 


? 


1, 2, 10 


18286-0959 




0.24 


200 


~15 


3.1 


5, 11 


18450-0148 


W43A, OH31.0+0.0 


0.8 


190 


50 


2.6 


6, 8, 13, 15,16, 20 


18460-0151 


OH 31.0-0.2 


0.11 


290 


~5 


6.8 


5, 11 


18596+0315 


OH 37.1-0.8 




59 






4 


19134+2131 




0.15 


100 


40 


8.0 


12, 14, 16 



fFull range of the line-of-sight velocities of maser emission. 
^Dynamical age of the jet. 

«fl: IBaud et al. 19791 2 : IBoboltz fc Marvel 20051 3: IClaussen et al. 200414: IDeacon et al. 20071 
5: |Deguchi et al. 2007[ 6: Diamond et al. 1985| 7: |Engels 2002) 8: IGenzel fc Downes 19771 
9: |Engels et al. 1986| 10: I Gomez et al. 19941 11: Imai et al. 2007 in p repa- 
ration, 12: |Imai, Sahai, fc Morris 2007[ 13: llmai et al. 20051 14: llmai et al. "20041 15] 
llmai et al. 20021 16: ILikkel et al. 19921 17: IMorris et al. 20031 18: ISahai et al. 19991 19: 
|Suarez, Gomez, fc Morata 2007[ 20: |Vlemmings et al. 2006] 

in the Galaxy, leading to estimating masses and birth points of the stars driving these 
water fountains. 

2. Lifetime and timing of the water fountain 

Table [T] gives physical parameters of the water fountain jets derived from single-dish 
and interferometric observations. Lengths of stellar jets traced by H 2 maser emission 
range over 6000 AU (in IRAS 16342-3814, IMorris et al. 20031 IClaussen et al. 2004p . 
which is still equal to or shorter than a radius of a circumstellar envelope of an AGB 
star (e.g-. IElitzur 1992]) . On the other hand, flow velocities of them (Vh ow ~ AVi os /2) are 
sometimes over 100 km s™ 1 . These give the dynamical ages of jets shorter than 100 years, 
only ^5 years in an extreme case (Imai et al. in preparation). In the time interval of 
100 years, a stellar jet traveling with a typical water fountain velocity of 100 km s _1 
reaches ^2000 AU from the central star, which corresponds to the typical size of an 
OH maser shell in a circumstellar envelope around an OH/IR star, at the final stage of 
energetic mass loss. The H 2 maser emission is excited near the tip of the stellar jet 
where high compression takes place in the shock at the interface between the jet and the 
ambient circumstellar gas. But because the ambient gas density declines with radius, the 
maser emission may be quenched at the distance where the post-shock gas density falls 
below a critical value, n H2 ~10 6 cm -3 (jElitzur et al. 19 89). 

Note that while optical lobes are observed in some water fountains (IRAS 16342—3814, 
and IRAS 19134+2131), the optical emission is not detected in others (W43A and 
OH12.8— 0.9), suggesting that the presence of water fountains may help in the identi- 
fication of new pre-PNe, which may be optically invisible due to heavy circumstellar ex- 
tinction. The optical lobes in IRAS 19134+2131 and IRAS 16342-3814 coexists with the 
H2O maser emission (|Sahai et al. 19991 IMorris et al. 20031 |Imai, Sahai, fc Morris 2007fl . 
It is expected that the lobes have been recently formed along the collimated fast jet, in 
which the circumstellar envelope around the evolved star has been penetrated by the jet 
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and starlight. These sources may be at the similar stage in the transition from an AGB 
to post-AGB star. 

Does the size of a circumstellar envelope determine the lifetime of a water fountain? 
We already know that H 2 masers can survive until photoionization of the circumstellar 
envelope starts to form a PN, such as K3— 35 ([Miranda et al. 200~lj) . In K3— 35, H 2 
maser emission is seen at the tips of the radio lobes as well as in the inner and denser 
part of the circumstellar envelope where where UV radiation from the star has not 
yet penetrated. However, the maser emission at the lobe tips of K3-35 is only weakly 
collimated, both morphologically and kinematically, implying that an initially collimated 
jet traceable by H2O maser emission is being destroyed. The dynamical age of the radio 
lobes is estimated to be 800 yr (jMiranda et al. 20"0T|) . which may correspond to an upper 
limit to the lifetime of a water fountain source but is one order of magnitude longer than 
the dynamical age estimated using VLBI observations of the maser spatio-kinematical 
structures. 

Here we speculate that the collimated jet may drill through the circumstellar envelope 
formed during the AGB phase and the tip of the jet may reach the outer surface of the 
envelope before the onset of photoionization of the envelope. The masers lying beyond 
the optical lobes, such as those in IRAS 16342—3814, might soon be quenched within a 
few decades or centuries owing to the decline of the gas density in the pre-shock gas as 
the jet propagates outward, which in turn leads to a loss of collimation and ultimately 
to a loss of the maser inversion as the post-shock gas density declines. 



3. Equatorial flows in the water fountain sources 

Highly collimated fast stellar jets as described above are expected to be driven by 
magneto-hydrodynammical process (e.g.. IBlackma net al. 200l| |Vlemmings et al. 2006[ ). 
However it is obscure whether the stellar object driving the jet is a single star or a 
binary system; one is an AGB star providing material to another through an accretion 
disk. Astrometry for H2O and 1612 MHz OH maser emission will provide an important 
clue to find each of the locations of the driving sources of the collimated jet and the AGB 
flow. Only W43A has SiO maser emission flNakashima fc Deguchi 2003[IImai et al. 2005| ). 
which provides a better opportunity of exactly find the location of the AGB star. 

Figure [1] compares the locations of the dynamical centers of the H2O maser jet and the 
SiO maser flow using a common extragalactic position reference. Although astrometric 
accuracy is much poor in the north-south direction for W43A located very close to the 
celestial equator, the dynamical centers are coincide within 4 mas (10 AU at 2.6 kpc). 
The SiO maser flow is well modeled by a biconically expanding flow with a low velocity 
(^20 km s _1 ) and a large opening angle flmai et al. 2005[) . These implies that both of 
the collimated jet and the biconical flow are driven in the region within 10 AU. 

Note that the spatio-kinematical structure of the SiO masers does not exhibit any 
rotation or contraction as suggested when mass accretion occurs from a disk. Figure [2] 
shows relative proper motions of the W43A H2O masers. We found a maser feature that 
is located close to the dynamical center and has a Doppler velocity close to the systemic 
velocity of W43A (Vl S r =34 km s^ 1 ) and a proper motion corresponding to a speed 
of ~30 km s -1 . The velocity vector is pointed to the different direction from the jet 
major axis. The low velocity maser feature is expected to be associated with an "equato- 
rial" flow, rather than an accretion disk. Such equatorial flows have recently been recog- 
nized in other water fountain sources, IRAS 18286-0959 and IRAS 18460-0151, as well 
( [Deguchi et al. 2007[ ). Figure [3] and H shows the spatio-kinematics of IRAS 18460-0151 
and IRAS 18286—0959, respectively. The former source has the fastest jet ever observed 
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RIGHT ASCENSION (J2000) 



Figure 1. H2O, SiO (v =1, J =1-0) maser emission, and 7 mm continuum emission in W43A. 
Different colors denote maser features' line-of-sight velocities, according to the color scale. Left 
Spatio-kinematics of the H2O maser features observed on 2002 April 3. A dotted line shows a 
pattern of a precessing jet appearing in the year 2002.3, which is modeled bv llmai et al. 20021 A 
symbol "X" indicates the dynamical center of the modeled jet. A thin dashed line indicates the 
axis direction of the H2O maser jet. Right Spatio-kinematics of the SiO maser features. Horizontal 
and vertical bars of individual maser features indicate uncertainties of feature positions in R.A. 
and decl. directions, respectively. A plus symbol indicates a roughly-estimated location of the 
dynamical center of the biconically-expanding flow found in the SiO maser kinematics. Roughly 
a flux density of a spot is inversely proportional to an error bar of the spot. 
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Figure 2. Relative proper motions of H2O masers in W43A measured during 2003-2004 with 
a mean proper motion subtracted. A cross indicates the location of the dynamical center of the 
precessing bipolar jet. Left the whole distribution of maser features. Right zoom-up view of the 
feature distribution in the area displayed in a dashed-line box in the left panel. 



with a speed over 300 km s . This source also has an equatorial flow with a velocity 
of ~30 km s _1 in the region close to the midpoint of the most blue-shifted and red- 
shifted maser features forming the high velocity jet. The latter source has many velocity 
maser components, only the most blue-shifted and red-shifted maser features are clearly 
separated each other and each of them is distributed in a small area. These imply the 
coexistence of a collimated fast flow and a slow outflow associated with either an AGB 
envelope or an equatorial flow driven by the fast flow as its sub-product. The existence 
of such equatorial flows may support a hypothesis that a single star can create both of 
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fast collimatcd and slow equatorial flows when the star is collapsing into a white dwarf 
(H. Washimi and S. Miyaji, in private communication). 



4. Ballistic corkscrew jets 

The morphology of water fountains may provide a unique clue to elucidating how to 
eject the material in the jets, leading to estimating the character of the central stars. Some 
of the water fountain jets exhibit spiral patterns. In particular, such a pattern is clearly 
seen in the distributions of H2O masers in W43A (|Imai et al. 20021 llmai et al. 2005]) . 
On the other hand, a "corkscrew jet" is found in optical emission in IRAS 16342—3814 
(jSahai et al. 2005)) . Recently we found twin spiral patterns of the H2O maser distribution 
in W43A, one of the spiral is seen in the axisymetrically opposite side of another (dotted 
line) with respect to the jet major axis (thin dashed line) (see figure [TJ However, curved 
motions of the maser features are not recognized. The right panel of figure [2] and figure [5] 
show proper motions of H 2 masers found in W43A and IRAS 19134+2131, respectively, 
in detail. In each source, we find an alignment of maser features whose direction has an 
offset from the major axis of the jet. On the other hand, proper motion vectors of maser 
features are still parallel to the jet axis. These suggest that the corkscrew jets found in 
these sources are formed by ballistic motions of maser features from the driving source of 
the jets, but whose directions are precessing with time. The alignments of maser features 
may be formed at the shock fronts interfaced with the ambient circumstellar envelopes. 



5. Location and motion of the water fountain in the Galaxy 

H2O masers associated with IRAS 19134+2131 (hereafter 119134) have measured 
their absolute coordinates with respect to the extragalactic position reference source 
J1925+2106 for a span over one year fllmai, Sahai, fc Morris "2007) ) . We therefore fit the 
temporal variation of the mean position of three maser features to a kinematical model 
consisting of an annual parallax and a constant velocity mean secular motion. Figure [6] 
shows the fitting result. We obtain a best-fit annual parallax, corresponding to a distance 
to 119134 D = 8.0^g'7 kpc. On the basis of the inferred heliocentric distance and the 
systemic secular maser motion, we estimate the location of 119134 in the Galaxy to be, 
{R,0,z) = (7Atoi kpc, 62°±5°, 0.65^;^ kpc, and the 3-D velocity vector in Galactic 
cylindrical coordinates to be, (V R , V e , V z ) = (3±g, 125^, 8±f|) [km s -1 ]. 

The estimated location is close to the Perseus spiral arm within ~1 kpc (c.f., Nakanishi & Sofuc 2006) 
Note that the estimated Galactic rotation velocity of 119134 (~125 km s" 1 ) is much slower 
than that implied by the Galactic rotation curve (~220 km s _1 , e.g., IClemens~ 1985, 
|Dehnen ~fe BinncyT998|. It implies that the I19134's orbit must have an orbit with large 
excursions in radius although 119134 is orbiting along a circle at present. 119134 is appar- 
ently orbiting as a member of the Galactic "thick" component rather than the "thin" com- 
ponent. Adopting the source velocity perpendicular to the Galactic plane V z =8 km s _1 , 
we find that the present Galactic altitude gives a lower limit to the travel time of 119134 
from the Galactic plane to the present position to be ~ 2.4 x 10 7 yr. This suggests that 
the stellar object of 119134 has a lifetime of at least this travel time and a stellar mass 
smaller than ~8.6 M (here we adopt a stellar lifetime to be oc M*/L* ex M~ 2 8 ) if 
the central star of 119134 was born in the Galactic plane as a member of the Galactic 
"thin disk" component. The statistical study of Galactic PNe bv IMancha do 2004 shows 
that the Galactic scale height of bipolar PNe such as 119134 is < z >=100 pc, suggesting 
that their progenitors should be high mass stars. Although the mass of I19134's progen- 
itor is therefore expected to be near the upper mass limit estimated above, the present 
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Figure 3. Distribution (filled circle) and rela- 
tive proper motions (arrows) of H2O masers in 
IRAS 18460-0151 during 2006 September-Novem- 
ber. Left Overview of the distribution. Right Zoom-up 
view of the area displayed with a dashed line box 
in the left panel. Different colors denote the fea- 
tures' line-of-sight velocities, according to the color 
scale above this panel. Relative proper motions of 
maser features are measured with respect to the 
maser feature located at the origin and biased in the 
north-west direction by ~30 km s _1 . 
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Figure 4. Same as figure [3] but for 
IRAS 18286-0959 observed on 2006 September 14. 
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Figure 5. H2O maser feature distribution in IRAS 19134+2131 (roots of arrows) and relative 
proper motion vectors of these masers (arrows) with respect to a position reference maser feature 
with a mean proper motion subtracted. 
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Figure 6. Mean motion of three H2O maser features in IRAS 19134+2131 with respect to the 
extragalactic position reference source J1925+2106. (a) Mean R.A. and decl. offsets on the sky 
with respect to the phase-tracking center, (b) and (c) Mean R.A. and decl. offsets against time. 
(d) Relative mean offsets with a mean proper motion subtracted. The observation code is shown 
alongside each measured position. 



galactic altitude of 119134 and its anomalously slow galactic rotation velocity are quite 
odd. These may imply an abnormal environment during the birth of I19134's central 
star. However, the estimated mass of 119134 still has a big uncertainty due to the poor 
accuracy of the vertical velocity, V z . 
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